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Abstract 

We present a microscopic, analytical theory describing a confined iV-electron 
gas in two dimensions subject to an external magnetic field. The number 
of electrons N and strength of the electron-electron interaction can be arbi- 
trarily large, and all Landau levels are included implicitly. For any value of 
the magnetic field B, the correlated iV-electron states are determined by the 
solution to a universal effective problem which resembles a ficticious particle 
moving in a multidimensional space, without a magnetic field, occupied by 
potential minima corresponding to the classical iV-electron equilibrium con- 
figurations. Introducing the requirement of total wavefunction antisymmetry 
selects out the allowed minimum energy iV-electron states. It is shown that 
low-energy minima can exist at filling factors v = 2 n+i where p and n are 
any positive integers. These filling factors correspond to the experimentally 
observed Fractional (FQHE) and Integer (IQHE) Quantum Hall effects. The 
energy gaps calculated analytically at v = | are found to be consistent with 
experimental data as a function of magnetic field, over a range of samples. 
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I. INTRODUCTION 

The problem of a highly-correlated, two-dimensional electron gas in an external magnetic 
field has attracted much attention in the past decade. Of particular interest is the micro- 
scopic origin of the observed fractions in the Fractional Quantum Hall Effect (FQHE)0H1. 
In the past few years, it has also been appreciated that many-body effects play a role in the 
formation of the gaps giving rise to the Integer Quantum Hall Effect (IQHE). As a com- 
plement to the experimental work in this subject, there have been many theoretical models 
proposed for both the FQHE and the IQHE. These range from field-theoretical treatments 
through to numerical, finite-size (N < 6) calculations. One of the most successful theoret- 
ical developments has been the proposal of trial wavefunctions by Laughlin and othersEIIHi 
to describe the interplay of wavefunction antisymmetry and electron-electron repulsion that 
effectively allows electrons in the lowest Landau level to form a highly correlated electron 
liquid. A related development by Jainl considers the construction of 'composite' fermions by 
attaching flux tubes to each electron - recent work on Chern-Simons field theories provides 
some support for such composite-fermion construction schemeJl@. The general problem of 
describing an iV-electron gas in an external magnetic field has recently taken on additional 
importance in semiconductor physics due to the fabrication of quantum dots containing a fi- 
nite number of electrons§~0. It is interesting to note that although the FQHE was originally 
observed in infinite two-dimensional electron gases (2DEG), it even persists in quantum dots 
containing a large but finite number N of electrons^. 

Given the fact that the underlying, microscopic iV-electron Hamiltonian is known, one 
could ask whether there exists an alternative, more direct way of understanding the nature 
of highly-correlated electron states without recourse to composite-fermion constructions, 
effective field theories, restrictions to lowest Landau levels and/or small numbers of electrons. 
The obvious stumbling blocks are that the electron-electron repulsion and the cyclotron 
energy are typically comparable in magnitude, and that iV-electron Schrodinger equations 
are generally intractable analytically. 

In this paper we pursue such an alternative approach starting with an iV-electron 



Schrodinger equation. We develop a microscopic, analytical theory describing correlated 
states of a confined iV-electron gas in two dimensions subject to an external magnetic field 
B. The number of electrons N and the strength of the electron-electron interaction can be 
arbitrarily large, and all Landau levels are included implicitly. We show that the description 
of iV-electron correlated states at finite B reduces to a universal effective problem which 
resembles a ficticious particle moving in a multidimensional space occupied by potential 
minima corresponding to the classical iV-electron equilibrium configurations. Introducing 
the requirement of iV-electron wavefunction antisymmetry selects out the allowed minimum 
energy iV-electron states. A possible connection with the FQHE and IQHE is then proposed. 
In particular, it is argued that low-energy minima can form at particular angular momenta 
corresponding to filling factors v = where p and n are any integer. These filling factors 
correspond to those observed experimentally in the FQHE and IQHE. 

The present theory suggests the following possible physical interpretation of FQHE and 
IQHE states. Consider an iV-electron wavefunction localized around a Wigner crystal (WX) 
configuration with total relative angular momentum J. At particular values of J, iV-electron 
wavefunctions localized around nearby defect configurations (i.e. WX plus defect which we 
shall denote as WXD) can co-exist; we note that the allowed values of J such that N- 
electron states can co-exist around WX and WXD simultaneously are severely restricted by 
the requirement of total wavefunction antisymmetry. At these common J values, which we 
shall denote as J = J m , hybridization of the iV-electron states centered on the WX and 
WXD minima can occur. This hybridization effectively allows the electrons in the WX solid 
to diffuse throughout the system via WXD defect states. The resulting delocalized 'liquid'- 
like iV-electron state has a lower zero-point energy - a gap therefore opens up between 
the liquid-like states at J = J m and other states at J ^ J m . For large iV the resulting 
liquid-like ground-states at J = J m have filling factors given by the well-known formulai 
v = N *2j~^ • ^ e fi n d ^ e v values at which these gaps arise are identical to those 
observed experimentally in the FQHE and IQHE. The energy gaps calculated analytically at 
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v = I are found to be consistent with experimental data obtained from a range of samples. 
Various other known features of FQHE states can also be reproduced. 

The model avoids discussion of one-electron properties such as Landau levels, and there- 
fore offers the possibility of a unified description of both the FQHE and IQHE based on 
a microscopic iV-electron Schrodinger equation. The formalism in this paper builds on an 
earlier model presented by us in Ref. 12. In particular, we conjectured in Ref. 12 that the 
classical minimum energy configurations play a crucial role in deciding the symmetry-allowed 
iV-electron correlated states in few-electron quantum dots. It was pointed out that the clas- 
sical minimum energy configurations for N < 6 all consist of N particles on a ring, while for 
N = Q additional minima occur Curiously, it is precisely at N = 6 that the magic num- 
ber J sequence of A J = N is broken. This idea was independently pursued by Maksym in a 
fascinating wajJll for N < 6 - the classical Eckardt frame was employed to study correlated 
few-electron dynamics and, in particular, the possible existence of 'liquid'-like states. We 
note that the term 'liquid'-like was introduced by Maksym to describe the loss of symmetry 
occuring when states corresponding to different classical minima are allowed to mix. This 
terminology will also be used in the present paper. We wish to emphasize that the model 
presented here is qualitatively different from an earlier theoretical approach of Kivelson et 
al.@ based on the so-called cooperative ring exchange. In short, we are suggesting here 
that FQHE states are the liquid-like states resulting from the hybridization of iV-electron 
wavefunctions localized around both crystal (WX) and crystal-plus-defect configurations 
(WXD). 

The outline of the paper is as follows. In Sec. II we present the microscopic iV-electron 
Schrodinger equation. The hyperangular coordinate system is introduced for the relative 
motion Hamiltonian. The problem then reduces to a 2N — 4 hyperangular equation (Sec. 
II. 1). In Sec. II. 2 the specific case of iV = 3 is outlined. This was discussed in detail in Ref. 
12 and is reviewed here since it is useful for visualizing the iV-electron results. In Sec. II. 3 a 
simplified hyperangular equation is obtained which is valid in the regime of strong electron- 
electron interactions and for large N. The characteristics of the lowest energy solutions are 
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discussed. Section III addresses the requirement of iV-electron wavefunction antisymmetry. 
Permutation symmetries of the iV-electron wavefunction become space-group operations 
in the multi-dimensional hyperangular configuration space. The states which will become 
ground states separated by a finite energy gap are found to correspond to filling fractions 
observed in the FQHE. Section IV obtains analytic estimates for the FQHE gaps at the 
fractions | as an example. These estimates are found to be consistent with experimental 
data over a range of samples, despite the fact that the results emerge from a simple one- 
dimensional, particle-in-a-box equation. Section V summarizes the results. 

II. MICROSCOPIC ^-ELECTRON HAMILTONIAN 

The analytic tractability of our model is made possible by a combination of a parabolic 
confinement potential and an inverse-square electron-electron repulsion potential. The 
parabolic confinement is known to be a reasonable approximation for many semiconduc- 
tor quantum dot samples^. For the case of a heterostructure (i.e. 2DEG) it mimics the 
effect of a positive background yielding an approximately uniform electron density in the 
large N limit (see Ref. 14). The (5/r 2 electron-electron interaction (fl > 0) is not unrealistic 
in quantum dots due to the presence of image charges in neighboring electrodes. Recent 
theoretical work suggestsEl that the true repulsive interaction between electrons in a quan- 
tum dot dot is more likely to be proportional to l/r n with n ~ 3 at large r and n ~ 1 at 
small r. In heterostructures, the electron-electron interaction is probably less affected by 
image-charge effects. However the general features of our results, which are based on the 
assumption that n = 2 for all r, should still be qualitatively useful. In particular, the occur- 
rence of FQHE in two-dimensional electron gases is not thought to depend crucially on the 
precise form of the electron-electron repulsion. Recent quantitative comparisonsQIHIl have 
indeed shown that the 1/r 2 and 1/r repulsive interactions yield iV-electron energy spectra 
with very similar features. Of particular relevance to the present theory is the finding that 
the classical minimum energy configurations for N electrons in a two-dimensional parabolic 
confinement potential seem to be almost identical for both 1/r and 1/r 2 interactions^. 
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The exact Schrodinger equation for N electrons with repulsive interaction (3/r 2 moving 
in a 2D parabolic potential subject to a magnetic field B (symmetric gauge) along the z-axis, 
is given by (if S p ace + ff sp m)# = E^; 

N (p. _ eAi)2 I Q 
^pace = E( \J + ) + E 

= V(j£- + -m*u 2 JB)\r t \ 2 + ^k) + y , 13 - (1) 

2 

where uJq(B) = + c<j c is the cyclotron frequency, and H spin = —g*^BBJ2i s i,z- The 
momentum and position of the z'th particle are given by 2D vectors and rj respectively; 
li is the z-component of the angular momentum. The exact eigenstates are written in terms 
of products of spatial and spin eigenstates obtained from H space and H spin respectively 
Eigenstates of H sp[n are just products of the spinors of the individual electrons and have 
energy E spin = g*fiBBS z , where S z is the ^-component of total spin and g* is the electron 
effective ^-factor. We employ standard Jacobi coordinates Xj (j — 1,2, ... , N) where Xi = 
jj J2j T j (center-of-mass coordinate) and X J>1 (relative coordinates) is given by 

X, = l^Hrj ~ J^r(ri + r 2 + . . . r^)] (2) 

together with the conjugate momenta Pj (see Fig. 1 for N = 3). The center-of-mass 
motion decouples, H sps , ce = if CM (Xi) + i/ re i({X j> i}), hence E spSLCe = E CM + E Ieh The 
exact eigenstates of Hqm and energies Eqm are well-knownl!. The non-trivial problem is 
to solve the relative motion equation H TC \tp = E Te ii/j. We transform the relative coordinates 
{Xj>i} to standard hyperspherical coordinates: X, = r(J\f = j +1 cosa^sinoje^ with r > 
and < ctj < | (a 2 = f). Physically, the hyperradius r is just the root-mean-square 
electron-electron separation. The exact eigenstates of if re i have the form ip Te \ = R(r)F(Q) 
where Q denotes the (2N — 3) hyperangular {#; a} variables; R(r) and F(Q) are solutions of 
the hyperradial and hyperangular equations respectively. The hyperradial equation is given 
by 

d 2 2N-3d 7 (7 + 2N-A) r 2 2m*(E rcl - hJf ) 1 

[ ^ + —r~ Tr 7 2 1 + j? ]R{r) = (3) 



where 1$ = Ti(m* loq(B))~ 1 and J is the total relative angular momentum. The parameter 
7 > and is related to the eigenvalue of the B and a; -independent hyperangular equation 
(see Sec. II. 1). Equation (3) can be solved exactly yielding 

E rcl = hcu ( B) [2n + 7 + N - 1] + (4) 

where n is any positive integer or zero and 

R(r) = [f VLl +N -*&)f% . (5) 

Equation (4) provides an exact and infinite set of relative breathing-mode excitations 
2hu (B)An for any N regardless of particle statistics and/or spin states. These quan- 
tum breathing modes were first reported in Ref. 12, and later confirmed by Geller et al.@ - 
the classical version of these modes for the Coulomb interaction was discussed in detail by 
Schweigert et al.S 

II. 1 Exact hyperangular equation for any N 

It remains to solve the B and a;o-independent hyperangular equation which is given by 

2m* B 

[G% + -^V(Q)]F(Q) = [7(7 + 2N - 4)]F(fl) (6) 

where 

_ d 2 [2N-6- (2N - A)cos2a N ] d 2 2 2 d 2 

Bv = -tt^-H — 7^ h sec aAr9 A r_ 1 - cosec a N — T . (7) 

oa N sinzajv oo>n o0 n 

The quantity V(ft) is given by r 2 Ya<j l r « — r jl -2 anc ^ om y depends on hyperangular coordi- 
nates Q = {8; a}. We emphasize that this hyperangular equation (Eq. (6)) is universal in 
that it is independent of the values of the magnetic field or confinement: solving Eq. (6) for 
7, and hence using Eq. (4), yields the complete solutions of the iV-electron Hamiltonian H 
for all magnetic fields and confinement strengths. Unfortunately the hyperangular equation 
does not admit exact solutions for 7. Sections II. 2-4 will consider various approximations 
to Eq. (6) which make the problem tractable. Because J remains a good quantum number, 
we can introduce a Jacobi transformation of the relative motion angles in particular 



and 

%i = i J —^m - + fis + ■ ■ ■ (9) 

where j = 3, 4 ... iV. We hence have one 9' variable, (N — 2) 0y] variables, (AT — 2) a 
variables and one hyperradius r giving a total of 2N — 2 variables as required for the relative 
motion. The exact eigenstates of H Te \ have the form ip = e tje ' R(r)G(Vl') where Q! denotes 
the (2N — 4) {^fji; ctj} variables excluding 9' . The term V(Q) is independent of 9' and will 
hence be written as V(fi'). It is useful to rewrite the eigenvalue of the hyperangular equation 
in terms of a new variable e as follows: 

e = |[ T ( 7 + 2N - 4) - J 2 - ^VTO] (10) 

where V(Q' Q ) is the value of V(ft') evaluated at the hyperangles corresponding to a particular 
classical, minimum-energy iV-electron configuration (Wigner molecule). Permuting electron 
indices will provide a set {f^} of symmetrically equivalent minima (SEM)@>0 with the same 
potential energy V(fi;) = V(Q f ) for all i (e.g. Q' Q and Q[ shown schematically in Fig. 2). 
Such SEMs have the same topological structure but cannot be transformed into each other 
by rot at tonsilS As will be shown in Sec. II. 2, there are two such SEMs for A^ = 3. The 
quantity e in Eq. (10) accounts for the contribution to the eigenvalue of the hyperangular 
equation without including either the contributions from the rigid-body rotational energy 
J 2 or the electrostatic potential energy ^r^-V(^l' ) of the classical minimum-energy config- 
uration. Physically therefore, e contains the zero point energy in f2'-space associated with 
the quantum-mechanical spread of G(Q') about the minima {f^}- The actual spread in 
G(Q') and hence e will depend on the total wavefunction antisymmetry requirement. This 
is illustrated for A^ = 3 in Sec. II. 2 and discussed for large A^ in Sees. II. 3, II. 4 and III. 
In general e > 0, e ~ (3^ where /i < 1, and e ~ J 5 where 5 < 2: these statements will be 
illustrated in Sec. II. 2 for A^ = 3. It is straightforward to show that the term 2 ^^ V(£l' ) 
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appearing in the definition of e is identical to [^|fj] 2 where Vd ass is the potential energy of 
the classical, minimum-energy iV-electron configuration, thereby recovering the expression 
given in Ref. 12. Note that Kiass oc P^uq(B) and that e (like 7) is independent of B and 
ujq. The exact relative energy for any iV can now be written as 

E rel = hw (B)[2n + ([N - 2] 2 + J 2 + -^V(Q' ) + + 1] + j_£ . (n) 

E ie \ only depends on particle statistics through e. As Ti — > 0, e — > and _E rc i — > Ki ass . 

The exact -E re i expression has an important consequence. The J-dependence of e is weaker 
than J 2 as J — > 0. Hence the term J%r in ^rei will dominate the J-dependence of _E re i for 
small J at fixed magnetic field u c . For states with J < 0, E Te \ will hence initially decrease as 
\J\ increases. On the other hand at large negative J, E rei will tend to h(uj (B) — ^)] J| and 
hence will increase linearly with |J| at a given u; c . This implies that E Tei has a minimum at 
a finite negative J for a given fixed magnetic field u c . This is the basic mechanism behind 
the tendency of an iV-electron gas to form ground states at increasingly large J values as 
the magnetic field is increased. As will be shown in Sec. II. 2 for N = 3 electrons, and in 
Sec. Ill for large N, only a subset of these J minima are permitted under the requirement of 
total wavefunction antisymmetry. These J-values are often called 'magic number' J-values 
in the context of few-electron quantum dots. In Sec. Ill we will show that the analogous 
'magic number' J states for a large- N electron gas constitute FQHE and IQHE states. We 
emphasize that so far our results are exact for any electron number N, electron-electron 
interaction strength /3, magnetic field uo c and parabolic confinement c^o- 

II. 2 Specific case of iV = 3 

This case was studied in Ref. 12. Here we will summarize the results since they are 
important for understanding the general iV case. For convenience we change variables from 
a, 9 (c.f. Fig. 1) to x,y where x = ln[tan(| — a)} and y = 6 — |. Since < a < |, 
hence —00 < x < 00 (N.B. —ti < y < n). We define p x = j-^ and p y = kj^. The exact 
hyperrangular equation (Eq. (6)) now takes the form 
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2 f , ftJcos(2tan e x )\2 



where 



(2 + cos(2tan 1 e x ) 3 9/ _-, „. 

Wa;, u; e) = m*/3h ; t — s —, 1 xxn 5 -snr(2tan x e ) 

v ' y ' y ^ L (cosec(2tan- 1 e a; ) + cot(tan- 1 e 2: )) 2 -3sin 2 ?/ 4 v 1 

+icos 2 (tan- 1 e s ) + ^-cos 2 (2tan- 1 e a; )] . (13) 
2 m*p 

Equation (12) represents the single-body Hamiltonian for a ficticious particle of energy e 
and unit mass, moving in the xy-plane in a non-linear (i.e. e-dependent) potential V(x, y; e), 
subject to a ficticious, non-uniform magnetic field in the ^-direction 

5 fic = ^[l-cos(4(tan-V))] . (14) 

.Bfic is independent of B and has a maximum of at x = for all y. For small x, 
Bfr c « — x 2 ). As x — > ±00, Bfi c — > 0. Note we have here chosen to highlight the 

Schrodinger-like form of Eq. (12); a simple rearrangement of Eq. (12) shows it to be 
hermitian with a weighting function sin 2 (2tan~ 1 e :c ). These results are exact so far. Figure 3 
shows the potential V(x, y; e) in the (x, y) plane. V(x, y; e) > everywhere. Minima occur at 
(0, 0) and (0, ±7r) where V(x, y; e) = (N.B. (0, tt) is equivalent to (0, — 7r)). Maxima occur 
at (ln\/3, ±|) in Fig. 3, where V(x, y\ e) — >■ 00. Since e > 0, these statements hold for any 
e. We now discuss the physical significance of these features. The classical configurations of 
minimum energy (Wigner molecule) correspond to the particles lying on a ring in the form 
of an equilateral triangle with Vd ass = u (B)[6m* /3]3 . There are two distinct configurations, 
i.e. two distinct symmetrically equivalent minima@, with clockwise orderings Q' Q = (132) 
and Q[ = (123) corresponding to (a, 9) = (?,±|). In (x,y) coordinates, these correspond 
to (0,0) and (0,7r) (equivalently, (0, — 7r)). Hence the classical configurations coincide with 
the minima in V(x, y; e) in Fig. 3 and the maximum in B^ c . As pointed out in Ref. 12, the 
formation of a Wigner molecule should therefore be favored by both large B^ c (i.e. large 
\J\) and deep V(x, y;e) minima (i.e. large /?, strong electron-electron interactions). 

Consider the limit of very strong electron-electron interactions (i.e. (3 — > 00). Since the 
tunnel barrier height between the two V(x,y;e) minima ~ f3, the ficticious particle sits at 
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one of these minima and the system is locked in one of the two classical configurations, e.g. 
Q' Q = (132) at (0,0). The tunnelling probability between the minima Q' Q and Q[ is zero. 
Tunnelling between the two minima implies a mixture of configuration (123) into (132) and 
hence interchange of the original electrons; in many-body language exchange effects arising 
from wavefunction antisymmetry are therefore negligible, e is small compared to m*f3 and 
Eq. (12) reduces to 



2 , 6m*/3 N i 11 T hu> c 



E iel = hu (B)[2n+(l + J 2 + —^)*+l] + J^ . (15) 

The energy E ve \ > Kiass since it includes the hyperradial zero-point energy (N.B. h — > 
yields E vcl -> Kiass and B &c 0). 

Next consider large but finite f3. The ficticious particle now moves in the vicinity of the 
minimum Q' Q (i.e. (x, y) ~ (0, 0)). The electrons in the Wigner solid are effectively vibrating 
around their classical positions. Expanding the potential V(x, y;e) about (0,0) to third 
order, the exact Eq. (12) becomes 



[ 2 2 2 x 2~ v 



+ ^ 2 7 + ~y x x 2 + -u 2 y y 2 ]G(x, y) = eG(x, y) (16) 



where uj 2 x = ( 3m A 13 + 2e) and uj 2 = 3 ™ 13 . This has the form of a single electron moving in 
an anisotropic parabolic potential, subject to a uniform magnetic field B^ c = Equation 
(16) is exactly solvable for e using a symmetric gaug (the energies are independent of the 
choice of gauge for -Bfi c ). As an illustration, we consider small e hence uj x ~ u y . The relative 
energy is then given by 

E rel = hu (B)(2n + [1 + J 2 + + 2(2n' + |/|' + 1)(J 2 + + 21' jf* + 1) 

. (17) 

The ficticious particle has its own set of Fock-Darwin (and hence Landau) levels^ labelled 
by n' and a ficticious angular momentum For large j3 and small n', I' and J, Eq. (17) 
yields an oscillator excitation spectrum with two characteristic frequencies \/2huJo(B) and 
2hu)o(B) representing shear and breathing modes of the Wigner molecule. 
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For smaller (3 (i.e. weaker interactions) and/or larger e (i.e. excited states), the tunneling 
probability between the V(x, y; e) minima Q' Q and Q[ in Fig. 3 becomes significant. The 
Wigner molecule begins to melt and wavefunction antisymmetry must be included. This is 
discussed further in Sec. III. As mentioned in Ref. 12, the resulting analytically obtained 
magic-number J transitions are found to be in good agreement with the numerical results for 
1/r interaction. We note that the analytic results become more accurate in the Wigner solid 
regime (e.g. large (3 or \ J\) while the numerical calculations become more computationally 
demanding. 

II. 3 Simplified hyperangular equations for arbitrarily large N 

For general N the hyperrangular equation (Eq. (6)) is (2N — 4)-dimensional. However 
in the Wigner solid regime (large (3 or |J|) the classical minimum energy configurations will 
still be important in determining e and hence E Ie \, just as for N = 3. Here we will consider 
the limit that the number of electrons is large (N >> 1). This is the limit of interest in 
the FQHE and in large quantum dots. Specifically, we will introduce in this section a series 
of approximations in order to simplify the exact hyperangular equation. At each stage, the 
corresponding simplified hyperangular equation is explicitly given. The resulting discussion 
is detailed - however we feel that this is necessary in order to justify the successively simpler 
(and more approximate) hyperangular equations. Each of these simplified hyperangular 
equations can be solved numerically: the complexity of the algorithms needed obviously 
decreases as more approximations are introduced. However, the goal in this paper is to 
obtain a simplified version of the hyperangular equation which can be treated analytically, 
but which is still based on a set of reasonable approximations. 

Figure 4 shows the classical ground-state configuration for iV = 230 electrons (black 
dots) in a parabolic quantum dot, as obtained by Bedanov and Peeters using a Monte Carlo 
algorithm!!. The rings are drawn as a guide to the eye. The number N of electrons is 
relatively small in the context of the N — > oo limit and hence the details of the ground-state 
configuration, particularly for larger rings, will be prone to edge effects. However, the inner 
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rings show a nearly hexagonal lattice as expected for the N — > oo limit. For the purposes of 
illustration we will therefore consider Fig. 4 as being representative of the N — > oo classical 
configuration. Consider the particular classical configuration fl' where the iV'th electron is 
near the center and the first electron is on the circumference of the droplet. This is shown 
schematically in Fig. 5. As discussed in Sec. II. 2 for N = 3, the fully quantum mechanical 
system will also lie near this configuration in f2'-space in the limit of very large f3. The j'th 
Jacobi coordinate is given by 

X^^l^.-R^) (18) 

where Rj-i = jri( r i + r 2 + . . . r?-i): this quantity Rj-i can be thought of as the 'center- 
of-mass' of the electrons 1, 2 . . .j — 1. For j >> 1, and for configurations as in Figs. 4 and 
5 where the electrons are evenly distributed around the origin, the quantity Rj-i will be 
small compared to the typical electron lattice spacing. In addition, the prefactor [^j-]^ — > 1 
for large j. Hence Xj — > r,- for large j. However, there is an exact identity for hyperangular 
Jacobi coordinates: X/^tv"] 2 = ^ Given that Xj — > as j — > N for large N, this implies 
that each term [^-] 2 << 1 for large j. Given the definition of the hyperangular coordinates 
stated earlier on, it follows that the hyperangles aj « 1 for j » 1. Hence to first order in 
Oij, we can approximate Xn = rsinaTv ~ rax- Similarly X^-i = rcosaArsina7v-i ~ ra^r-i 
and, more generally, Xj m raj. To summarize, for configurations similar to that shown in 
Fig. 5, we have the approximate result X^- ks raje 10 ^ in the N » 1 and j >> 1 limit. 
There are two points to note: although we need both j » 1 and N >> 1, j can still be 
an order of magnitude less than N. Second, the error introduced by assuming sinaj <x,- 
is still reasonably small even for j = 2 (recall a.2 = f ~ 1-57 as compared to sina 2 = 1- 
To remain consistent within our approximation, we will take a 2 — 1 instead of | in what 
follows) . 

This approximate form for Xj leads to an interesting simplification of the exact hyper- 
angular equation. The small-angle (i.e. ctj « 1) limit of Eq. (6) yields: 
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N ^2 n 2 

E ~2W V ? + ^WW) = 2^ 7(7 + 2N ~ 4)F(fi) (19) 

J =2 

where 

V 2 - — — — — (20) 

J 0(\'j acj ih\ , (\j i)Oj 

is the two-dimensional Laplacian for a ficticious particle with position (aij,0j) in polar co- 
ordinates, the potential energy term 

v(q) = v(n') ~ Yl \ a i ei6j - a r ei6j ' I" 2 ( 21 ) 

i<i' 

and -F(fi) = e lje ' G(Vl'). This equation is a good approximation for j -+ N with TV >> 1 
but becomes worse as j — > and/or AT — > 0. (Recall a 2 = 1 hence the sum can start 
from j = 2 as shown). However this is sufficient for the purposes of this paper since we are 
interested in states that evolve within the bulk of the iV-electron droplet as opposed to those 
at the edge. Physically, Eq. (19) describes a set of N — 1 ficticious particles moving on a 
two-dimensional plane subject to a two-body inverse- square interaction, in the absence of a 
magnetic field. It is interesting to note this transformation of having replaced an iV-particle 
problem in a magnetic field with an N — 1 particle problem without a magnetic field seems 
reminiscent of composite fermion constructions at half-integer filling fractions. The effective 
Schrodinger equation in Eq. (19) carries the following constraint: the exact hyperangular 
identity Ef =2 [^} 2 = 1 implies £ 

a j ~ 1- This may complicate any attempt at a solution 
using a 'renormalization' type of approach, such as the setting up of a recursion equation 
relating 7 for iV particles to 7 for N — 1. 

It is more useful to view Eq. (19) in the context of a single ficticious particle moving in a 
multi-dimensional space containing potential miminima corresponding to the various classi- 
cal minimum energy configurations {f^}. This directly connects the iV-electron problem to 
the N = 3 problem discussed in Sec II. 2. As discussed in Sec. II. 1, a Jacobi transformation 
can be undertaken on the {9j} variables. In particular, 

0[i] = [y^fM-e,--i) (22) 
14 



where j = 3, 4...JV and 0j_i = ^2(^2 + #3 + . The quantity 0j_i represents 

the average of the angles 9j where j = 2,3. . .j — 1. For j » 1 the quantity Qj-i will 
be approximately a constant, 0, since the j — 1 particles are evenly distributed about 
the origin in a given f2- configuration (recall Figs. 4 and 5). In addition, the prefactor 
[jrf] 5 — *■ 1 f° r large j. Hence 9^ — > — for large j, neglecting terms of order (-^). With 
F(Q) = e iJ0, G(W), Eq. (19) further reduces to 

e + + 4 [ 4 + ^i ,2) + wn ' )IG(sr) = 

n 2 

" 7(7 + 2N — 4)G{Vl') (23) 



2m 



Again this equation is a good approximation for j — > iV but becomes worse as j — > 0. The 



{9[j];aj} manifold carries the following constraints: ~ 1 and Xig^ ~ 0. The latter 



d _ d 



»L 

condition is an approximate identity for large iV and is obtained by combining Y ^ (J 
(this is an exact property of any Jacobi transformation) and -J^- ~ + j^Zi^r- This 
new condition hence reflects the fact that the total relative angular momentum is only 
associated with the 9' variable: there is no additional contribution to the relative angular 
momentum contained within the Q' dynamics. These approximate constraints allow us to 
make a further simplification of the hyperangular equation as follows. Using the approximate 
identity J2 a ] ~ 1 we can define an average hyperangle a ~ N~^. We will therefore replace 
the term Y ^ (j^pp m Eq. (23) by Y (n-i) 2 ~ ^ assuming large N. We can hence 
rewrite Eq. (23) in the form 

tg-^ ( S + ^ + k [ i + ^t4 d + ;3|l ' (sy) - ,/W)l|G(sy) = 

7) 2 9m*R 

^[7(7 + 2iV - 4) - J 2 - -^y(fi' )]G(fi') . (24) 

Although its derivation has involved approximations, Eq. (24) merits some discussion since 
it elucidates several of the statements made in Sec. II. 1. The right-hand side is just -4re. 
Given that Y a | ~ 1 an d Y ~ 0, the J-dependence of e will tend to be weaker than J 2 
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as claimed earlier. Note that since ~ 1, the moment of inertia / of a given classical 
configuration in fi-space, treated as a rigid body, is just m* . Hence the classical rigid-body 
rotational energy ~ which is precisely the term appearing in the right-hand side 
of Eq. (24). This then justifies the statement made in Sec. II. 1 that e excludes the classical 
rigid-body rotational energy. The term V{Q! Q ) denotes V(Q') evaluated at a given classical 
SEM equilibrium configuration Q f = Q' . We emphasize that V(Q' ) = V(f2-), i.e. same 
potential energy for all SEMs. Since il' is a minimum, the difference term |V(f2') — V(fio)] 
can be expanded around fl' . The leading terms will be quadratic in 9y] — 9y ] etc. Hence e 
does indeed describe the zero-point energy associated with the spread in G(Q') around the 
classical minima, as claimed in Sec. II. 1 and shown explicitly for N = 3 in Sec. II. 2. This 
point is further discussed below for large N. 

The hyperangular equation Eq. (24) is now simpler, however it is still not quite in a form 
which makes it amenable to analytic calculation. This final step can be achieved with the 
following considerations. Given the two approximate constraints a ? ~ 1 an d J2 tS~ ~ 0> 
the term involving should be small as compared to the term involving \-Mr, and 

hence will be neglected. Furthermore just as for N = 3, we are initially considering the 
quantum mechanical solution near a given classical minimum Q' , i.e. (3 — > oo. Hence the 
term -^i-^p- can be approximated by -^r-§^r where aj is the value of otj at Q' = Q' . The 
fact that SIq is a minimum suggests that the leading order expansion of [V(fi') — V(fio)] 
will involve terms like (0^ — 9[jo]) 2 and (ctj — «jo) 2 f° r & h j but not cross terms: this was 
demonstrated explicitly for N = 3 earlier where V(x, y) was found to be a function of x 2 
and y 2 but not xy. This implies the following simplification for Q' ~ Q' Q : [V(Q') — V(Q' Q )] ~ 
J2j[v(aj — ctjo) + w (@[j] ~ Q\jo])] where both v and w have a minimum at Q' = Q' , i.e. at 
ctj = ajo and 9yj = 9\jo]. The approximate separability of the potential suggests that the 
hyperangular function G(Q') can now be written as f({ctj — otjo})g({0\j] — 0\jo]}) where the 
functions / and g are peaked around Q! = Q' . This was shown to be true explicitly for 
N = 3 where / and g turned out to be gaussians (harmonic oscillator wavefunctions). Since 
Qq is still a minimum point for large N, f and g will retain their gaussian-like character for 
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general N. We will therefore write g({9\j] — %o]}) ~ Tlj 9j(9\j] ~ 9\jo]) where gj(9y] — 9[j ]) is 
a function peaked around the minimum coordinate 9^ = 9y j. The hyperangular equation 
is now fully separable into an equation involving {aij} 

N h 2 d 2 Id 
^ ~2m^^do^ + '^'da^ + ~~ a j°^f({ a j ~ aj0 ^ = E <*f({ a i ~ a i°}"> ( 25 ) 

together with the following equations for each 9yy. 

h 2 d 2 

' ~ 2m* a 2 d9 2 + ^ Wj ^ ~ 6 m)\9j{9\j] ~ d [j0]) = e J g j (9 U] - 9 m ) . (26) 

3 ^ \J ] 

The relation between e, E a and ej is as follows: = Eg + E a where Eg = J2j e j- The full 
expression for the relative energy hence becomes 

1m* B 1m* i flu: 

E ml = huj (B)[2n+([N-2] 2 + J 2 + ^V(%) + ^[Ee + E a }^ + l] + J^ . (27) 

Since Eqs. (25) and (26) have a Schrodinger-like form with E a and ej as eigenvalues re- 
spectively, we will refer to these two quantities as 'energies' even though this is not strictly 
correct terminology. 

II. 4 Characteristics of the low-energy solutions G(Q') 

So far we have considered the solutions near a particular minimum configuration Q' , i.e. 
we have considered very large (3 just as we did initially for iV = 3. Very large j3 implies 
that G{0!) will be peaked around one of the SEMs, e.g. around Q' Q . In the limit of zero 
tunnelling between SEMs, the solution G(Q') ~ f({ctj — ®jo})g({9\j] ~ %o]}) centered around 

will be degenerate with the identically localized solutions centered at all other SEMs {f^}. 
These localized functions can be thought of as atomic-like orbitals in Q' space. In particular 
there will be a set of orbitals associated with each SEM f^. The corresponding coordinates 
and hyperangular equations describing these solutions are identical to those obtained earlier 
in Sec. II. 3: however the spatial ordering of the electrons for the various Q\ minima will 
necessarily change; for example electron N will not necessarily be close to the center. Using 
the usual variational argument for Schrodinger-like equations, the lowest energy (i.e. lowest 
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E a and Eg) solutions of Eqs. (26) and (27) will be those with the minimum number of 
nodes. 

For large but finite (3 there will be a small but finite tunnelling between the various 
minima {f^}, hence the complete solution G(Q!) will be more correctly described as a linear 
combination of the atomic-like solutions, just as in a single-particle tight-binding model. 
Furthermore for iV > 6, as noted earlier, there will be additional classical minima which are 
not topologically equivalent; again borrowing from the language of molecular physicUJi 
these minima are termed symmetrically inequivalent minima (SIM). These SIMs are local 
minima in f2'-space which are often just slightly higher in energy than the SEMs {fli}. In 
the large N limit, these minima correspond to defect states in a hexagonal crystal. Fisher, 
Halperin and Morfil showed that a Wigner crystal with a localized defect (WXD) can be 
quite close in energy to the perfect Wigner crystal (WX). This finding was recently verified 
in the context of N electrons in a two-dimensional parabolic quantum dot by BoltonEl and 
Bedanov et al.@. These authors all found that the global minimum for the classical N- 
electron system tends towards a hexagonal crystal as iV — > oo, as expected for the Wigner 
crystal (WX). However configurations corresponding to a Wigner crystal with single defects 
(WXD) are only slightly higher in energy. In the language of the present paper the WX 
represents the SEMs while the WXD represents the SIMs. Although the SIMs are not true 
global minima, the complete solution G(Q') should certainly include finite mixing with them. 
This is particularly true since the 'nearest-neighbors' of a given SEM in fi'-space are SIMs. 
This is simply a consequence of the fact that translation between two adjacent SEMs in 
f2'-space requires interchange of at least two electrons, while translation between a given 
SEM and its nearest SIMs requires only slight electron distortion. Each SEM minimum Qi 
will have p defect states as its nearest neighbors in f2'-space - we denote these nearby SIM 
minima as {£^. a } where a = 1, 2 . . .p (c.f. Fig. 2). 

The resulting wavefunction G(Q') will therefore resemble a tight-binding LCMO (Linear 
Combination of Molecular Orbitals) wavefunction where each 'molecule' consists of 'atomic' 
orbitals on one of the SEM minima f2j mixed with 'atomic' orbitals on each of its nearby 



SIM minima. We emphasize that N = 3 has no SIM minima. iV = 6 is the smallest N 
having SIMs. The SIMs for N = Q consist of a six-member ring configuration while the 
SEM's contain a five-member ring plus one electron at the center00. For general N, the 
low-energy solutions should therefore be reasonably well-described by 

SEM SIM 

G(fi')~ E E C ^/(K-^;a})^({%-%;a]}) • (28) 
i a 

It is well-known from elementary tight-binding theory that the lowest-energy states are 
'bonding' wavefunctions of s-orbitals. In the present context, we similarly expect the lowest 
energy G(Q') to have as few nodes as possible (i.e. it will be gaussian-like around each of 
the SEM {fl'j} thereby resembling an s-orbital); it will also correspond to the coefficients 
Ci- a being identical for each i (i.e. it will resemble a 'bonding' state). 

III. FERMION STATISTICS, MAGIC NUMBERS AND FILLING FRAC- 
TIONS 

So far we have not introduced the requirement that the total iV-electron wavefunction 
be antisymmetric. In this section, we will show that it is precisely this requirement that 
produces the observed FQHE filling factors for large N. 

It is useful to first discuss the effect of antisymmetry in the case of N = 3 electrons before 
considering large N. For three spin-polarized electrons, if) must be antisymmetric under 
particle interchange i <->• j. The hyperradial part R(r) is invariant; particle permutation 
operations in (r!,r 2 ,r 3 ) become straightforward space-group operations in the (x,y) plane. 
For small (x,y), 1 <-> 2 is equivalent to (x, y) — > (x, y + n) with 9' — > 9' + |; 1 <-> 3 is 
equivalent to (x, y) — > (x,y — n) with 9' — > 6' + | ((x,y) represents (x, y) rotated by 4^); 
2 <-> 3 is equivalent to (x, y) — > (x, y + ir) with 9' — ► 9' — | {(x,y) represents (x, y) rotated 
by — -f-)- The solutions G(x, y) of Eq. (12) with the lowest possible e and hence lowest E TC \ 
at a given u c , should be nodeless in the vicinity of (0, 0) (c.f. ground state in the parabolic 
potential with n' = = /' in Eq. (17)). However the above symmetry requirements forbid 

2 J 

such a nodeless solution unless e t7T ~ = 1. Therefore the only symmetry- allowed solutions 
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G(x, y) which are nodeless are those where J is a multiple of three, as observed in numerical 
calculations for N = 3 electrons with a 1/r interaction. It is important to note that this 

2 J 

condition, i.e. e ln ~ = 1, just arises from combining the effect of any two sets of particle 
interchanges i j. For N = 3, two sets of particle interchanges correspond to rotations 
of a given SEM: this can be seen simply as follows. Consider a given SEM in Fig. 3, e.g. 
Q' Q = (132). Interchanging 1 <-> 2 and 2 <-> 3 yields the same SEM, i.e. (132), rotated 
anticlockwise by Hence combinations of two sets of particle interchanges merely rotate 
the Wigner molecule without involving a transformation from one SEM to another, i.e. 
without moving from (132) to (123). Hence in order to obtain the 'magic' angular momentum 
values for N = 3, it is sufficient to consider the subset of particle interchanges from the S3 
permutation group which correspond to point-group rotations C3, i.e. those which do not 
correspond to translations between SEMs. This result is discussed by Maksym in Ref. 15 
following earlier work on molecules by Wilson@. Maksym also argued for N = 3 that the 
remaining permutations which correspond to translations between SEMs, and hence the 
effect of tunnelling between SEMs, represents a small perturbation which does not affect the 
magic J- values. In contrast, for N = 6 where topologically distinct classical configurations 
coexist, the tunnelling between SEMs and SIMs plays a crucial role in determining the J 
values of the low-energy ground states of the system. In particular, Maksym pointed out 
that tunnelling between SEMs and SIMs should be most favorable when the SEM and SIM 
configurations have a common J value. This is consistent with analogous ideas in single- 
particle tight-binding theory, where the overlap matrix element (and hence bandwidth) is 
larger between s orbitals than between s and p orbitals. Maksym conjectured that the 
resulting tunnelling might lead to 'liquid'-like states with a lower overall energy. 

These considerations motivate us to follow a similar strategy for N electrons. In particu- 
lar, we will show that considering just a subset of particle interchanges of Sn corresponding 
to rotations of rings within the Wigner molecule (WX) and Wigner molecule plus defect 
(WXD) is sufficient to determine the magic J- values corresponding to the observed FQHE 
filling factors. As for N = 3, we focus on the vicinity of a given SEM, e.g. Q' Q . Following 
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the discussion in Sec. II. 3, interchanging <-> r 3 - is relatively straightforward for i, j >> 1 
since X^ raje 1 ^ ~ r^. Neglecting terms of order (-^), it just corresponds to ctj <-> ctj 
and 0j <-> 0j. The derivation of the transformation rules including terms of order (-^), is 
straightforward but tedious. Just as for iV = 3, however, it turns out in what follows that 
we do not need to consider individual i <-> j transformations. 

III.l Spin-polarized system 

Consider the classical configuration Q' Q shown in Figs. 4 and 5. For large j the electrons 
can be thought of as forming an approximately ring-like structure. Counting the number 
of rings from the center outwards, the first ring contains 6 electrons, the second contains 
approximately 12 and so on. We first focus on a 'typical' ring without any defects: it will 
contain a large, even number N m of electrons (approximately 6m electrons where 6m » 1) 
but these electrons will have an index j » 1, i.e. we are not considering rings near the edge 
of the iV-electron droplet. We are going to consider just the subset of all particle interchanges 
i «-> j which are equivalent to rotations of this m'th ring. Since Xj ~ raje l9j ~ r 3 -, all 
members of the ring have approximately the same a, i.e. {aij} = a. m for all j in ring m. 
Hence interchanging two members of the ring just involves a transformation between their 
9j coordinates. Since all members of the ring have a similar environment and the same aij, 
the potential energy Wj and hence gj in Eq. (26) will have the same form for all j in the 
ring m. As for N = 3, the lowest-energy solutions should be those with gj nodeless; gj 
will be an approximately gaussian function of 9\j] centered around 9y ]. Hence we can write 
gj = g m for all j in ring m. We now rotate the electrons in the ring, and hence the ring itself, 
by an angle ■^ L . Since Xj ~ rj this corresponds to 9j — > 9j + jf-. The g m functions are 
nodeless and (in a given ring) identical, hence the product Ylj9j{%] ~ %o]) for j i n ring m 
can be replaced by Y[j 9m(6\j] — %o])- The transformation keeps the system within the subset 
of all SEMs corresponding to the same ring ordering, i.e. just as for iV = 3 the rotation 
operation in real space becomes a space-group operation in fi'-space which translates the 
system between SEMs. Recall that G(Q') for minimum-energy states should resemble a 



21 



'bonding' linear-combination of s-like orbitals (i.e. an approximately gaussian ^-dependence 
around the various SEMs (Eq. (28))). The coefficients Cj ;a in the expression for G(Q') in 
Eq. (28) will therefore be identical for all the SEM minima corresponding to this same 
ring ordering. Since G(Q') corresponds to a linear combination of identical orbitals with the 
same coefficient, the overall effect of the transformation on G(Q') due to ring rotation will be 
quite small; we simply move between this subset of SEMs {^-} each of which has the same 
local orbitals. In contrast, the effect on the 9' variable is relatively important; it follows from 
Eq. (8) that 9j -> 9j + corresponds to 9' -> 9' + A9' where A9' = jV m ^-|^ = 
The total function F(Q) = e ije 'G{Q') hence becomes e iJAe 'F(Q). Given that N m is an 
even number, rotation of the m'th ring by j^- necessarily corresponds to an odd number of 
interchanges i j. If we assume the electrons are spin polarized, the spatial part of the 
iV-electron wavefunction must be totally antisymmetric and hence the overall phase change 
must equal e l7r ( 2n+1 ) where n is any integer. Denoting the J value as Jwx we therefore obtain 
the condition 

Jwx = \(N-l)(2n+l) . (29) 

Importantly, this criterion for Jwx is independent of m and hence holds for all rings m. In 
other words, this criterion guarantees that the iV-electron wavefunction has the correct per- 
mutational symmetry under the subset of all permutations of N electrons which correspond 
to ring rotations. Note that Jwx must be an integer. 

Now consider the Wigner crystal plus defect (WXD). Fisher, Halperin and Morfil showed 
that the lowest energy defect states correspond to interstitial defects, i.e. an extra electron 
sits on an interstitial site in the otherwise perfect crystal. Single vacancies have a higher 
energy. Following Fisher et a!., there are two types of interstitial site, 'centered' and 'edge' 
interstitials, and these are by far the most predominant type of defect at finite temperatures. 
In our model, these defects can be created by introducing an (iV + l)'th electron which forms 
the defect. There are several reasons why this is reasonable. First, the alternative scheme 
of allowing one of the N existing electrons to form the defect would create an interstitial 
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plus vacancy; following Fisher et al. the total energy of such a defect is approximately three 
times larger than a single interstitial. Second, creation of such an interstitial-vacancy pair 
would involve a transformation of both 9 and a coordinates within the iV-electron f2-space. 
Third, the definition of the 9\j] variables (see Eq. (9)) is independent of the coordinates of 
electron N + 1. Hence the iV-electron Q' coordinate system is essentially unchanged by the 
presence of the extra electron. The hyperangular function G(Q') for the N electron system 
can therefore be compared directly to the corresponding hyperangular function for the N + l 
electron system when projected onto the iV-electron Q' space. We wish to consider the effect 
of this defect on ring m. Following Fisher et al.@ the distortion of the crystal will be well- 
localized around the defect. In terms of the hyperangular coordinates, part of the local 
crystal distortion will be subsumed in the coordinate r and the effect on the hyperangles 
a and 9 of the electrons in ring m will be relatively small unless the defect lies in ring m. 
Assume the defect lies in ring m = ntd- The antisymmetry condition obtained above for the 
perfect crystal (Jwx) will still be approximately valid for all rings with m ^ ma- In ring 
m d , there are now an odd number of electrons N m + 1. The rotation 9j — > 9j + N 2n +1 now 
corresponds to 9' -> 9' + A9' where A9 f = (N m + 1) (JVm + 1)JV = § (N.B. we now have an 
(A^ + 1) electron system). Because of the odd-member ring, rotation corresponds to an even 
number of % <-> j interchanges. The hyperangular function for the N+l electron system 
of crystal plus defect, when projected onto the original iV-electron fi'-space, is essentially 
unchanged - the original N electrons are only slightly distorted by the presence of the defect. 
Hence the overall phase change e lJAe ' must equal e l2wn ' where n' is any integer. Denoting 
the J value as Jwxd we therefore obtain the condition 

Jwxd = Nri . (30) 

Again, this criterion for Jwxd is independent of m and hence holds for a single defect located 
in any ring m. Also, Jwxd must be an integer. These two criteria, taken together, therefore 
guarantee that the A-electron wavefunction has the correct permutational symmetry under 
the subset of all permutations of N electrons which correspond to ring rotations, both for 
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the perfect crystal (WX) and the crystal plus defect (WXD). For the perfect crystal, we can 
consider N to be odd since each ring contains an even number of electrons, plus there is one 
electron at the center. Combining the two conditions for Jwx and Jwxd we hence see that 
the WX and WXD have the following common J values: 

J m = l -N{N -\){2n + l) (31) 

where n is any integer. Converting these J m -values into filling factors using the formula 
v = ^77—^, which is valid for large N, yields 



2 J n 

1 



(32) 



2n+ 1 

This coincides with the principal series of FQHE fractions, i.e. | and |. The value v — 1 
will be discussed below. As an illustration we consider the case of iV = 201 electrons. The 
allowed Jwx values are 100 x 1, 100 x 3, 100 x 5 etc. while the allowed Jwxd values are 
201 x 1, 201 x 2, 201 x 3 etc. It is clear that common J values are given by J m = 100 x 201 x 1, 
100 x 201 x 3 etc. and hence v — \, \. 

III. 2 Spin-unpolarized system 

We have so far generated the J m values for a system of spin-polarized particles. Next we 
consider the opposite limit of a spin-unpolarized system, i.e. N + = iV_ where N = N + + AL 
and N must therefore be an even number. The arguments will be more approximate in 
this case, but we believe will still contain the essential physics. Consider a 'typical' ring as 
before. Let this ring, m, contain N m electrons where N m » 1: the ring will typically have 
4p- up-spins and ^f- down-spins. Due to the Pauli principal keeping like spins apart, we will 
assume that on the average the ordering corresponds to the alternating sequence up-spin- 
down-spin repeated around the ring (See Fig. 6(a)). Rotation of the ring to a topologically 
identical configuration now involves a rotation of all the electrons in the ring by an angle 
(iV/2) ' * ,e ' we nave to r °t a te through twice j^. The rotation 9j — > 9j + j^jj) corresponds 
to 9' — > 6' + AO' where Ad' = N m (JVm/ 2)V-i) = W^) ' Since Nm is an even number > can 
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either be odd or even. Rotation of the m'th ring by therefore corresponds to either an 
even or odd number of interchanges % <-> j for both spin-up and spin-down electrons. Hence 
the total number of interchanges of like spins is always even. The overall phase change must 
therefore equal e l27vn where n is any integer. Denoting the J value as J W x we therefore 
obtain the condition 

Jwx = \{N-l)n. (33) 

Again this criterion for J W x is independent of m and hence holds for all rings m. Now 
consider the Wigner crystal plus defect (WXD) with the defect in ring m. The defect 
corresponds to an extra electron which can either be spin-up or spin-down. There are now 
an odd number of electrons N m + 1. We now have, for large JV m , that AO' rs 2(N m + 
1) , N 2 +^ N = 77 for the N + 1 electron system. Because of the odd-member ring, rotation 
now corresponds to an overall odd number of i «-> j interchanges. This is because either the 
spin-up interchanges are odd while the spin-down are even, or vice versa. Hence the overall 
phase change e lJAe ' must equal e l7T ( 2n '+ 1 ) where n' is any integer. Denoting the J value as 
Jwxd we therefore obtain the condition 

Jwxd = \n(u' + i) . (34) 

Again, this criterion for Jwxd is independent of m. These two criteria, taken together, 
therefore guarantee that the spin-unpolarized iV-electron wavefunction has the correct per- 
mutational symmetry under the subset of all permutations of N electrons which correspond 
to ring rotations, both for the perfect crystal (WX) and the crystal plus defect (WXD). 
Combining the two conditions for integer values of Jwx and Jwxd we hence see that the 
WX and WXD have the following common J values: 

J m = ^N(N -l)(2n + l) (35) 



where n is any integer. Converting these J m values into filling factors yields 

2 



2n+ 1 
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(36) 



This coincides with the second series of FQHE fractions, i.e. f , |, f etc. and hence suggests 
that the ground states at these fractions will be spin-unpolarized in the absence of Zeeman 
energy, in agreement with earlier finite-size numerical calculations (see p. 63 of Ref. 3). 
Interestingly the series also reproduces the IQHE value v — 2. As an illustration we consider 
the case of iV = 200 electrons. The allowed J W x values are 199 x 1, 199 x 2, 199 x 3 etc. 
while the allowed Jwxd values are 50 x 1, 50 x 3, 50 x 5 etc. Common J values are given 
by J m = 199 x 50 x 1, 199 x 50 x 3 etc. and hence v — |, | etc. 

III. 3 Partially spin-polarized system 

Consider a partially spin-polarized system. First we will take i.e. N + = 3AL where 
N = N + + N_ and iV is again even. Consider a 'typical' ring as before. Let this ring, m, 
contain N m electrons where N m » 1: the ring will typically have up-spins and ^ 
down-spins. Due to the Pauli principal keeping like spins apart, we will now assume that on 
the average the ordering corresponds to the sequence up-spin-up-spin-up-spin-down-spin 
repeated around the ring (see Fig. 6(c)). Rotation of the ring to a topologically identical 
configuration now involves a rotation of all the electrons in the ring by an angle jj^j^, i.e. 
we have to rotate through four times The rotation 9,-, — > 9j + /r ,? 7T ,., corresponds to 
9' -> 6'+A6' where A9' = N m jj^j^-^ = ^^j. Since N m is an even number, ^ and 
are either both odd or both even. Rotation of the m'th ring by therefore corresponds 
to either an even or odd number of interchanges i <-> j for both spin-up and spin-down 
electrons. Hence the total number of interchanges of like spins is always even. The overall 
phase change must therefore equal e' l27rn where n is any integer. Denoting the J value as 
Jwx we therefore obtain the condition 

Jwx = \{N-l)n. (37) 

Again this criterion for Jwx is independent of m and hence holds for all rings m. Now 
consider the Wigner crystal plus defect (WXD) with the defect in ring m. The defect can 
either be spin- up or spin-down. There are now an odd number of electrons N m + 1. We 
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now have, for large N m , that AO' ps 4(N m + 1) ( N J+i) N = for the Af + 1 electron system. 
Because of the odd-member ring, rotation now corresponds to an overall odd number of i <-> j 
interchanges. This is because either the spin-up interchanges are odd while the spin-down 
are even, or vice versa. Hence the overall phase change e tJAe ' must equal e m ( 2n ' +i ) wnere n > 
is any integer. Denoting the J value as Jwxd we therefore obtain the condition 

Jwxd = \N{n' + i) . (38) 

Again, this criterion for Jwxd is independent of m. These two criteria, taken together, 
therefore guarantee that the partially spin-polarized A^-electron wavefunction has the cor- 
rect permutational symmetry under the subset of all permutations of A" electrons which 
correspond to ring rotations, both for the perfect crystal (WX) and the crystal plus defect 
(WXD). Combining the two conditions for integer values of Jwx and Jwxd we hence see 
that the WX and WXD have the following common J values: 

J m = l -N{N -l)(2n + l) (39) 

where n is any integer. Converting these J m -values into filling factors yields 

4 



v = 



2n+l 

4 4 



(40) 



This coincides with the fourth series of FQHE fractions, i.e. |, | etc., hence the theory sug- 
gests that the corresponding ground states at these fractions will be partially spin-polarized 
in the ratio of spin-up to spin-down of 3 : 1 in the absence of Zeeman energy. Again this is 
in agreement with finite-size numerical calculations (see p. 63 of Ref. 3). As an illustration 
we again consider the case of A^ = 200 electrons (N.B. we have 4 as a factor since the ratio 
of spin-up to spin-down is 3 : 1). The allowed J W x values are 199 x 1, 199 x 2, 199 x 3 etc. 
while the allowed Jwxd values are 25 x 1, 25 x 3, 25 x 5 etc. Common J values are given 
by J m = 199 x 25 x 1, 199 x 25 x 3 etc. and hence v — |, | etc. 

Next we take i.e. A+ = 2AL where A^ = A^_ + AL and A^ is again even. Consider a 
'typical' ring as before. Let this ring, m, contain N m electrons where N m » 1: the ring 
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will typically have up-spins and ^ down-spins. Due to the Pauli principal keeping like 
spins apart, we will assume that on the average the sequence corresponds to up-spin-up-spin- 
down-spin repeated around the ring (see Fig. 6(b)). Rotation of the ring to a topologically 
identical configuration now involves a rotation of all the electrons in the ring by an angle 
jjTj^j , i-e. we have to rotate through three times jf-. The rotation 9j — > 6j + jw^m 
corresponds to 9' — > 9' + A9' where A0' = N m ^ N ^ N _^ = (jv-i) • Since N rn is an even 
number, 2 -^ n - and ^ are either both odd or both even. Rotation of the m'th ring by j?- 
therefore corresponds to either an even or odd number of interchanges i <-> j for both spin-up 
and spin-down electrons. Hence the total number of interchanges of like spins is always even. 
The overall phase change must therefore equal e t2nn where n is any integer. Denoting the J 
value as Jwx we therefore obtain the condition 

Jwx = \{N-l)n. (41) 

Again this criterion for Jwx is independent of m and hence holds for all rings m. Now 
consider the Wigner crystal plus defect (WXD) with the defect in ring m. The defect can 
either be spin- up or spin-down. There are now an odd number of electrons N m + 1. We 
now have, for large N m , that A9' ~ 3(A^ m + 1) ^ 2 ^ N = for the N + 1 electron system. 
Because of the odd-member ring, rotation now corresponds to an overall odd number of i <-> j 
interchanges. This is because either the spin-up interchanges are odd while the spin-down 
are even, or vice versa. Hence the overall phase change e tJA9 ' must equal e' 1 ' 2 "'" 1 " 1 ' where n' 
is any integer. Denoting the J value as Jwxd we therefore obtain the condition 

Jwxd = \N{n' + I) . (42) 

Again, this criterion for Jwxd in independent of m. These two criteria, taken together, 
therefore guarantee that the partially spin-polarized iV-electron wavefunction has the cor- 
rect permutational symmetry under the subset of all permutations of iV electrons which 
correspond to ring rotations, both for the perfect crystal (WX) and the crystal plus defect 
(WXD). Combining the two conditions for integer values of Jwx and Jwxd we hence see 
that the WX and WXD have the following common J values: 
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J m = l -N{N -\){2n + l) (43) 



where n is any integer. Converting these J m values into filling factors yields 

3 



2n+ 1 

3 3 



(44) 



This coincides with the third series of FQHE fractions, i.e. |, | etc. and predicts the 
corresponding ground states to be partially spin-polarized in the ratio of spin-up to spin- 
down of 2 : 1 (or vice versa) in the absence of Zeeman energy. We note that there is 
an alternative system that also yields the filling factor series v = ^pi- m particular, this 
fraction emerges from considering a fully spin-polarized system but now considering rotation 
of two rings simultaneously. These two states with iV + : AL = 2:1 and N- = respectively, 
probably compete to become the ground-state depending on the value of the Zeeman energy 
(and hence magnetic field). Interestingly, finite-size studies have shown that the v = | state 
is indeed partially polarized for B < 15T in the ratio N + : AL = 2:1 but fully polarized 
for B > 15T (see p. 160 of Ref. 3). 

The above arguments can be extended straightforwardly to consider N + — {p — 1)AL 
where N = N + + AL and p is any integer. In this case, the corresponding filling fraction 
becomes v = . We now focus on the filling factor v = 1 (IQHE). We have already 
shown that this state emerges from the two following series: v = for n = in a fully 
spin-polarized system, or v = -^-^ for n = 1 in a partially spin-polarized system. In fact 
for a given spin polarization iV + — {p — 1)AL with p odd, the factor u = 1 will always 
arise, i.e. by choosing n = e y l . It is reasonable to expect these states to compete to 
become the ground-state. This is consistent with recent findings that a gap exists at v — 1 
(IQHE) even in the absence of Zeeman splitting. The possible coexistence of a manifold of 
partially spin-polarized states is also consistent with the idea of macroscopic spin-textures 
near v = 1; in particular, taking a linear combination of partially spin-polarized states 
enables the construction of localized 'wave-packets' of spin - we conjecture that the resulting 
spin-textures may be related to skyrmions. 
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We note that the above filling factors emerged from requiring that the WX and WXD 
had a common angular momentum value. The defect was considered to be an interstitial 
electron. It turns out that, as far as the common J m values are concerned, we could also have 
considered the defect to be a vacancy. The product (N — 1)(N — 2) would have appeared 
throughout this section instead of N(N — 1); in the large N limit, both products yield 
~ N 2 and hence the same filling factors. Hence at J = J m , the WX can coexist with 
a WXD where the defect is either an interstitial electron or a vacancy. As noted earlier, 
however, such vacancies do have a higher energy and hence are less likely to occur at finite 
temperatures. 

IV. ANALYTIC CALCULATION OF FQHE ENERGY GAPS 

In this section we will give an analytic calculation for the energy gaps associated with 
the FQHE and IQHE states, i.e. for J = J m . The calculation is approximate since it relies 
on the various approximations made in Sec. II. 3. However our goal is to find whether the 
gaps predicted by our model are in fact consistent with the observed FQHE gaps and also 
to identify trends in the energy gaps with filling fraction, magnetic field etc. 

In Sec. II. 3 we obtained an approximate expression for the relative energy E m \ (see 
Eq. (27)). This energy depends on E a and E e . In Sec. Ill we argued that the important 
criterion characterising the magic-number J-values was that the crystal (WX) and defect 
(WXD) can both have the same J value given by J = J m ; this leads to a large derealization 
energy due to increased WX-WXD tunnelling in f2'-space. In the language of single-particle 
tight-binding theory, the resulting energy gap between states with J = J m and J ^ J m 
arises from the hybridization of the G{0!) solutions peaked around, for example, Q' and 
Q' Q . a at J = J m . This hybridization hence yields a low E Te \ because of the corresponding 
derealization of G(Q') at J = J m , i.e. a reduction in zero-point energy. Here we will obtain 
an analytic expression for this energy using a simple model for the effect of derealization 
and show that the resulting gaps are consistent with experimental findings. 

As pointed out earlier, a state of a given negative J will have an energy minimum at 
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a finite magnetic field u c . As u c increases, the value of J at which the energy E re i(J) 
has a minimum will increase. We will calculate the energy difference between a state with 
J = J m and competing low-energy states with J given by J m ± = J rn ±5 at a given cu c , where 
5 « J m . The lowest energy state with J = J m is given by Eq. (27) with n = 0: 

£ rcl (J m ) = huo(B)[([N - 2f + J 2 m + -j^V{%) + -j^[E 6 {J m ) + E a {J m )]Y + 1] - J m ^ 



(45) 



while that with J = J m ± is given by: 



Er Cl (J m ±) = hw (B)[([N - 2] 2 + [J m ± 5] 2 + 

2 ^-V{%) + ^[E e (J m± ) + E a {J m± )\)^ + 1] - [J m ± 5}^ . (46) 

As discussed in Sec. II. 1, the angular momentum is negative for low-lying energy states, 
hence we have included the minus sign directly into these expressions, i.e. J m and J m ± 
are positive numbers. Consider the state with J = J rn . This state can coexist as both a 
crystal (WX) and a crystal plus defect (WXD), in contrast to the state J m ±. The transition 
WX-^WXD will involve a distortion of the ^'-coordinates. This distortion in f2'-space 
corresponds to a spreading or 'derealization' of the function G{Q!) along these directions. 
There is therefore a reduction in 'localization' energy going from a state with J ^ J m (i.e. 
WX and WXD cannot coexist) to a state with J = J m (i.e. WX and WXD can coexist). 
We will hence write: 

E e ( J m ) + E a ( J m ) + (5V{%) = E (47) 

while 

E e (J m± ) + E a (J m± ) + (3V(n' ) = E + A (48) 

where A represents the increased localization energy of state J m ± as compared to J m . Note 
that the potential energy minimum /3V(fl' Q ) is a constant term throughout. In the Appendix 
we show that typically A ~ iV 2 while Eq ~ A^ 3 . Hence in the limit of iV >> 1, we have 
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Eq » A. We also recall from Sec. Ill that J m ~ A^ 2 . It follows by expanding out Eqs. 
(45) and (46) in the limit iV >> 1 that 



AE„ = E re i(J m ±) - E re i(J m ) « huj {B) 



=F^. (49) 



We are interested in the large N limit since our goal is to calculate the FQHE gaps: hence 
we will choose the confinement ujq << uj c which yields uq(B) ~ i y. Let us first consider 
the filling factor v — | hence J m = 3 jV ( j ^~ 1 ) _ Substituting into Eq. (49) we obtain the 
approximate expression for the energy gap 

1 m*A 

AEi ^huj c . (50) 

There are several points to note about this expression for the energy gap at v = |. 

(i) Given that A ~ A^ 2 for large N, the expression is independent of the electron number 
N. It is also independent of the strength of the parabolic potential u>o- For a given uj c 
corresponding to the filling factor v — |, we can therefore take the thermodynamic limit 
iV — > 00 and yet still maintain a fixed average electron density by choosing appropriately 
small values of c^o- Our expression for the energy gap at filling factor v — g which was derived 
in terms of J values for a fixed- N system therefore also holds for an infinite two-dimensional 
electron gas of fixed density. 

(ii) The expression for the energy gap does not exhibit a direct dependence on the value 
of the electron-electron interaction {3. We do emphasize, however, that throughout most of 
this paper we have assumed that (3 is large enough to be able to neglect tunnelling between 
SEMs. Hence we can only conclude that the absolute value of (3 does not directly affect 
the energy gap AE for sufficiently large (3. This is consistent with experimental findings 
that the gap can be remarkably sample-independent^. Below we will mention how a weak 
dependence on (3 will arise if one considers smaller values of (3. 

(iii) The energy gap appears to be approximately linear 4n the magnetic field. Most previous 
theoretical studies conclude that the dependence resembles 55. As we will show below, the 
linear dependence is in reasonable agreement with experimental data, particularly at lower 
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fields. However, we will later discuss how a weaker, non-linear dependence, i.e. B x where 
x < 1, can eventually arise at larger B in our model. 

(iv) The energy gap does not depend on 5 to first order in (j-)- This independence of 5 is 
important since it implies that the energy gap exists between state J m and all other states 
with J in the vicinity of J m . We know from the discussion at the end of Sec. II. 1 that, over 
a given range of magnetic field, the states which compete to become the ground state will 
be those of similar J. Hence we expect the energy gap arising to exist over a small but finite 
range of magnetic field, as observed experimentally. 

(v) The expression for the gap AE can be made applicable to situations where the lowest- 
energy excitation involves spin-flips, by adding AE spin where AE spin is the difference in total 
spin energy between the excited state and the v — | spin-polarized ground state. However, 
we are interested in the lowest-energy excitations, and hence will take AE spin = since 
AE S pi n > for the fully spin-polarized initial state at v — |. For other fractions, the term 
A will have an indirect dependence on the spin configuration since, for a mixed spin system, 
the fl' space is really coupled to the spinor space by the antisymmetry condition, hence 
G(Q') is actually a two-dimensional vector. A detailed discussion of spin-reversed excited 
states will be given elsewhere. 

(vi) Given the discussion in (v), together with the fact that A has a weak but finite de- 
pendence on J rn (see later), we can conclude the value of A will generally be different for 
different fractions. Gaps at other fractions in the v — | series are discussed below. 

We will now attempt to derive an approximate analytic expression for A. Consider 
the iV-electron system at J = J m . As discussed in Sec. Ill, the system can exist as 
both a crystal (WX) and a crystal plus defect (WXD), in contrast to the state J m ±. We 
will argue in what follows that the transition WX— >WXD involves a significant fractional 
distortion of ^-coordinates. This distortion along the {9y]}-axes in f2'-space corresponds to 
a spreading or 'derealization' of the function G(Q') along these directions, thereby giving 
rise to a finite A. In Sees. II and III, we argued that the low-energy SIMs near a given 
SEM consisted of a single electron defect placed at an interstitial site within the hexagonal 
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crystal. As noted at the end of Sec. Ill the defect could also be a vacancy, although the 
corresponding SIM would have a higher energy. Consider a defect placed in ring to which 
contains N rn particles with N m » 1, and let the defect be sited between j and j + 1 
(c.f. Fig. 5). Classically, the system moves to a nearby SIM in f2'-space (c.f. Fig. 2). In 
particular the defect will cause a distortion of the coordinates of particle j. In order to 
calculate the maximum possible distortion (and hence derealization available as a result 
of the hybridization between the SEM and SIM) we will consider the particular SIM in 
which only particle j moves to accomodate the defect. In principle, both the aj and Oy] 
coordinates will be modified, thereby 'sharing' the effect of the distortion. However, with the 
defect placed between j and j + 1 in ring to, the distortion of particle j will mainly be along 
the By] direction. The idea that the most important effect of the defect is the distortion of 
the ^-coordinates, is consistent with the following considerations. Consider any particle j' 
nearby to the defect with coordinates 9[y] and ay. Let the nearby defect cause a distortion 
of a in all directions. Hence the new coordinates of the particle j' are approximately ay + a 
and 6\~n + — . The relative distortion caused by the nearby defect is hence ~ — while 
~ a .,g ■ For N >> j', ay changes slowly with j' and is unchanged if f lies in ring 
to. However 9^ ranges from to 2n within ring to and hence the fractional change A0y] 
can be significant. 

We will therefore consider A as arising as a result of the difference in E e for J m as 
compared to J m ±. The loss of 'localization' energy of the state J m as compared to J m ± can 
be significant along the {0ui}-directions, i.e. Eg can differ appreciably depending on whether 
the function G(Q') is localized (i.e. J = J m ±) or not (J = J m ). Consider Eq. (26) for the 
'energy' ej associated with a particle j in ring to where J ^ J m (i.e. WX and WXD cannot 
coexist). The equation resembles a one-dimensional Schrodinger equation in 8\y for a mass 
aj moving in a potential Wj(9[y —8\jo\)- The function Wj will have a minimum at 9y] = 9y j, 
whereas there are maxima at 9[y = 9\jo] ± jf- as particle j approaches particle j ± 1. For 
N = 3, recall Fig. 3 where moving along the y-axis (#-axis) at fixed a = a (i.e. x — 0) 
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produced a minimum at y = and maxima at y = ±|. We can approximate ej using a 
simple one-dimensional particle-in-a-box model: we will assume that Wj is a flat-bottomed 
potential with infinite walls at 9y] = 9y ] ± j^-. The width of the box is therefore given by 
dj = The energy ef^ is hence approximately given by 

j 2m*a] [47r] 2 ' 1 J 

Note that the width of the wavefunction gj, and hence the localization of G(Q') along 9y], 
is characterized by aj. This expression (Eq. (51)) for e^p implicitly assumes that the 
electron-electron interaction j3 is large; for smaller values of /3, the particle-in-a-box energy 
should pick up a weak dependence on (5. Throughout this paper, however, we will use the 
approximate analytic form given in Eq. (51). At J = J m there is distortion of G(Q') since 
WX and WXD can coexist. The defect can occupy any interstitial site in the crystal; each 
defect position produces a distinct SIM. There exists a SIM (e.g. f2J, -) in which the defect is 
placed next to particle j, say between particles j and j + 1 as used earlier. We can see that 
there is one such SIM associated with each 6y\ coordinate. If, as discussed above, we consider 
the main distortion as occuring on the 9y] coordinate of particle j, then the associated SIM 
lies on the 9y] axis. In this case the coexistence of the SEM £1' and the SIM Q'q.j causes an 
increase in the effective box-width, aj — * aj + 5a j. The energy along the 9y] direction is now 
given by 



O (J) jr 

(Jm) (J) . ° e j x (J) M n da j-i 

e) ~ e) + -Q^- Sa i ~ e ) t 1 " 2-^] • (52) 



Although there are many such SIMs associated with each SEM, we know that the defects 
have a low density at the temperatures of interest. The system also requires a large time 
to tunnel between these SIMs since each SIM describes a different defect position in the 
crystal: diffusion of the defect between sites will be slow at low temperatures. It is reasonable 
therefore to suppose that each SEM hybridizes with just one of these SIMs at any time. The 
average loss of localization energy of state J m as a result of distortion due to a nearby SIM 
is therefore obtained by averaging over all (N — 2) 9y] coordinates. Hence using Eq. (52) 
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where x represents an average of the quantity x over all (N — 2) of the 9u\ coordinates. 

We could also try to obtain expressions for AE at other fractions. Consider v = |. 
Equation (50) now has the factor | replaced by |. Assuming as a crude approximation that 
the values of A are the same, we obtain the result that AEi : AEi ~ 0.6 for samples at a 
given magnetic field. The literature tends to put this ratio at about 0.3 — 0. ME. One could 
also try to evaluate AE for other fractions where the ground-state J m and/or the excited 
states J m± 

are thought to have spin-reversed electrons. Such a calculation needs a more 
careful estimation of A, as discussed earlier. Here we will just provide a rough estimate 
by considering, as before, excitations which do not change the total spin component (i.e. 
AE spin = 0). We will choose v — | in order to compare with the results for v — |. There are 
two effects of considering v — | instead of v — | when re-deriving the expressions obtained 
in this section. The prefactor in Equation (50) increases by a factor of 2, while ^ decreases 
by a factor of 2. This decrease in — arises as follows. Recall that like spins in the spin 
unpolarized case are separated by twice the angle of the spin-polarized case. The 'unit-cell' 
size in fi'-space is determined by the separation between neighboring SEMs, and is therefore 
twice as large for the unpolarized case. The effective box size must therefore also be twice 
as large, i.e. the effective a* value is now 4r- instead of 4r-. Hence the final result is that 

J iVm iV m 

AE2 ~ AEi . Although we do not attach too much importance to this result because of the 
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complications of spin, it is interesting to note that experimentally the gaps for v — | and 
v — I are found to be similariSil as will be shown below. We note that allowing for spin- 
flip excitations at v — | may reduce the overall gap at low B since AE sp i n can be negative if 
the ground state is not fully spin-polarized. Hence the total gap may be negligible for u — | 
at low B. This feature is also seen experimentally. The same argument concerning effective 
box size should also be approximately true for the other partially spin-polarized fractions in 
the I series. Consider v = |. The prefactor in Equation (50) increases by a factor p, while 
the effective box-size also increases by the same factor. The net effect is that the expression 
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for the gap is similar for all fractions | where p = 1,2, 4,5 etc. Hence the energy gaps for 
the | fractions measured across a range of samples should all fall on approximately the same 
curve as a function of magnetic field. 

Equation (53) presented an analytic expression for A which, to the level of approximation 
employed, did not depend on the ground state J m value. Examination of the more accurate 
versions of the hyperangular equation presented in Sec. II. 3, suggests that A should actually 
have a weak but finite dependence on J m (recall, for example, Eq. (23) or (24)). In the 
N = 3 study in Sec. II. 2, we found that increasing J did indeed increase the localization of 
G(Q') around the classical minima. For large N, as discussed in Sec. II. 3, the dependence 
on J m should be weaker; however the resulting hybridization between a given SEM (i.e. 
WX) and nearest-neighbor SIMs (i.e. WXD) should also decrease slightly as J m increases. 
Decreasing hybridization will reduce the value of A as J m increases. It is interesting to 
analyze the effect of this reduction in A for the separate situations of (a) a given sample 
over a range of v values, and (b) different samples at a given fixed v. In case (a), the number 
of electrons N is fixed. As the value of J m increases, the value of u c at which this J m value 
represents the ground state must also increase. If A decreases as J m increases, then A must 
also decrease with increasing u c . This reduction in A as uo c increases, if sufficiently large, will 
make AE — > 0: we suggest that this may be related to the predicted formation of a Wigner 
solid (i.e. gapless excitations) at very high magnetic fields. In case (b), fixed v means that 
increasing J m requires an increase in N (recall that J m = JV ^~ ) ■ If A decreases as J m 
increases, then A must also decrease as N increases. However, a fixed value of v means 
that increasing N requires an increase in magnetic field. Hence A decreases as u c increases. 
As mentioned earlier, this will tend to weaken the linear magnetic field-dependence of the 
theoretical gaps at fixed v to lo x c with < x < 1. Such a sub-linear dependence is consistent 
with recent experimental data at high fields^. In the estimates of the gap discussed below, 
however, we take as a first approximation the form of A presented in Eq. (53). Consequently 
the calculated gaps for a given v always increase linearly with magnetic field. 
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We now proceed to discuss appropriate values of [^] and hence calculate the gaps. The 
precise value for [^] will depend on the details of the crystal plus defect system (WXD). 
Here we will suggest reasonable lower and upper estimates, and argue that the particular 
value to be used will depend on the degree of disorder in the experimental samples. In 
particular, we will argue that the lower estimate is appropriate for disordered (i.e. lower 
mobility) samples while the upper estimate is appropriate for pure (i.e. higher mobility) 
samples. 

First we consider the lower estimate for [— ]. Fisher et al. found that the maximum 
distortion for a vacancy defect was about 12%, but that the value for the interstitial defect 
was 'considerably larger'EH. If the sample contains a significant impurity concentration, 
it is likely that interstitial electrons will have difficulty in diffusing through the iV-electron 
system. The function G(Q') will therefore have a restricted derealization for kinetic reasons. 
In the absence of interstitial defects, the derealization would be determined solely by the 
vacancies. We will therefore take the value of 12% as a lower estimated bound for [— ]. From 
the Appendix we have that e ~ V^^X and hence A ~ 0.4^-4- Substituting this into the 
expression for the energy gap, we obtain AEi ~ 0.014huj c meV and hence AEi ~ 0.lGhu} c 
Kelvin. Given that lmeV= 1.728£[T] (Tesla) for GaAs, we obtain AEi ~ 0.275[T] Kelvin. 
Hence at B = 20T, AEi ~ 5.5 Kelvin. Figure 7 compares this lower estimated bound of the 
energy gap at v — | (and hence v — | as explained above) with early experimental results 
obtained by Boebinger et al.0H over a range of relatively impure samples (i.e. significant 
impurity concentration). The agreement is surprisingly good; however we emphasize that 
our calculation is obviously fairly crude. Apart from improving the expression for AEi given 
in Eq. (50), one could do a better job in calculating the localization energy A from Eq. (53). 
Such improvements would almost certainly render the calculation of energy gaps within the 
present model numerical, although such calculations would be more straighforward than the 
original alternative of an iV-electron diagonalization. Various numerical improvements will 
be presented in a future publication; the goal of the present paper is to pursue a purely 
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analytical theory. 

Now we turn to an upper estimate for [— ]. We assume that the sample is pure and hence 
there is no kinetic reason for ignoring interstitial defects. We first recall our physical picture 
of the topology of the interstitial defect, stated earlier in this section. We consider a ring 
m containing N m electrons, including particles j — 1, j and j + 1 (c.f. Fig. 5). The angle 
between particle j and j + 1 is ■^ L , and the angle between particle j — 1 and j + 1 is jf-. 
We let the defect lie in the ring between particles j and j + 1. We are looking for an upper 
estimate on the value of [— ], hence we will assume that the only particle which moves to 
accomodate the extra electron is particle j. As before, the angle between j — 1 and j + 1 
is still but now there are two particles (j and the defect) within this angular range. In 
an equilibrium state (i.e. a SIM) the three angles between j — 1 and j, j and the defect, 
and the defect and j + 1, are all equal to t^-- The distortion of the effective box size for 
particle j will be determined by the angle between j and j + 1. Hence an estimate of the 
average distortion [^] is ^-(| — 1) divided by jf-, which gives |. Our upper estimated 
value of [&] is therefore 0.33. Fi gure 8 compares both this upper bound and the lower 
bound obtained earlier to experimental data obtained by the Oxford and AT&T groups for 
a range of relatively pure, high mobility samplesil@. The experimental values lie between 
the two bounds. This consistency between the present theoretical results and experiment 
lends support to our interpretation of the effect of sample purity. 

V. CONCLUSIONS 

A microscopic theory describing a confined iV-electron gas in two dimensions subject to 
an external magnetic field, was presented. The number of electrons iV and strength of the 
electron-electron interaction can be arbitrarily large. For any value of the magnetic field B, 
the correlated iV-electron states were shown to be determined by the solution to a universal 
effective problem: this problem resembles that of a ficticious particle moving in a multidi- 
mensional space, without a magnetic field, occupied by potential minima corresponding to 
the classical iV-electron equilibrium configurations. 
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A possible connection with the Fractional (FQHE) and Integer (IQHE) Quantum Hall 
effects was subsequently proposed. In particular, it was shown that low-energy minima 
can arise in the large N limit at filling factors v = where p and n are any positive 
integers. The energy gaps calculated analytically at v — | were found to be consistent with 
experimental data as a function of magnetic field, over a range of samples. Various other 
known features of FQHE and IQHE states were found to emerge from the present theory. 

While it is obviously extremely difficult to calculate many-particle energy gaps etc. ac- 
curately using an analytic approach, we hope that the general qualitative trends and orders 
of magnitude provided by the model will be useful in understanding the fascinating but 
complex field of highly-correlated iV-electron systems. We also hope that the model may 
begin to shed some light on the connection between the two limits of few-electron correlated 
states in quantum dots, and the infinite two-dimensional electron gas. 
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APPENDIX 

First we derive an approximate expression for e and hence A. The exact hyperangular 
identity Ef= 2 [v] 2 = 1 im P lies £ 

Oij ~ 1. The moment of inertia / of the (N — l)-particle 
system in fi-space, treated as a rigid body, is therefore approximately just m*. For large 
N, the density of particles will be approximately uniform: the moment of inertia for such 
a uniform disk is just ^(N — l)m*R 2 ~ ^Nm*R 2 , where R represents the disk radius and 
(N — l)m* is the total mass. Hence R 2 ~ ^, the average density of particles is ^~ 2 ^ ~ 
and the average particle-particle spacing is ~ Now consider the sum over energies ef^ 

from Eq. (51): 

Ve (J> ~y h 2 n 2 [N m } 2 

The quantity otjo for particle j in ring m is approximately m times the average particle- 

i 

particle separation: i.e. ctjo ~ m&^-. Replacing the sum over j by a sum over the rings m, 
and using the approximate result that there are 6m particles in ring m, yields 

> e) ' ~ > 6m - ~ > m 55 

Y m 64m*vrm 2 8m*vr Y 

where m max is the maximum ring number. m max is given approximately by the disk radius R 
divided by the particle-particle separation. Hence m 2 max ~ ^. For large m max , Y^ max m ~ 



Wmax- Hence 



2>f~™^. (56) 
^ 16m*7r 2 

We require the average value, e, with the average taken over all j. There are N — 2 such 



j coordinates, hence for large N we have e ~ 2 l n f 2 as claimed in Sec. IV. Given that 



2 JU2 



16m* 7r 2 



A ~ 2[— ]e (57) 



we obtain 



; Ja^h 2 N 2 
A ~ — 58 

1 a J 8m*7r 2 
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and hence A ~ N 2 as claimed. 

Second, we investigate the general iV-dependence of E . Following Eqs. (47) and (48), 
we assume that E is dominated by the classical potential energy at the SEM Q' = Q' , i.e. 
E ~ f3V(Vt' ). This is consistent with our assumption throughout the paper of considering 
configurations close to the classical minima. Hence 

E °~'%\^'°>"-«r^' <59) 

Replacing the denominator by a 2 yields E ~ f3= J2j<j> 1 an d hence E ~ f3N^{N — 1)(N — 
2). For large N, E ~ (3^- and hence E ~ A^ 3 as claimed. We note that while this 
derivation is crude, the final expression for E Q is not actually used in the calculation of 
energy gaps. The only result used is the conclusion that E » A for large N. 
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Figure Captions 



Figure 1: Jacobi coordinates for the N = 3 electron system. Reading clockwise, the classical 
configuration for three repulsive particles Q' Q = (132) corresponds to (a, 9) = (f , f ) (i.e. 
(x, y) = (0,0)); the other classical configuration Q[ = (123) corresponds to (a, 9) = (|, — |) 
(i.e. (x,y) = (0, 7r) or, equivalently, (0, — ir)). 

Figure 2: Schematic diagram showing a portion of the hyperangular fi'-space. Two sym- 
metrically equivalent minima (SEMs - labelled as f2-) are indicated by larger circles, while 
nearby symmetrically inequivalent minima (SIMs - labelled as f2^. a ) are indicated by smaller 
circles. 

Figure 3: Contour plot of ficticious potential V(x, y;e) in the (x, y) plane for the N = 3 
electron problem. The two symmetrically equivalent minima (SEMs - Q f and Q[ from 
Fig. 1) are shown. Minima in V(x, y;e) occur at (0,0) and (0, ±n) (i.e. at the classical 
configurations). Maxima occur at (lnv^, ±f ), where V(x, y; e) — > oo (i.e. electrons 2 and 3 
or 1 and 3 coincident). V(x, y; e) is positive and finite everywhere else. The same qualitative 
features appear for all e = 5 is used as an illustration). 

Figure 4: Ground-state classical configuration for 230 electrons (black dots) calculated nu- 
merically using a Monte Carlo simulation. Figure adapted from Fig. 2 of Bedanov and 
Peeters (Ref. 14). Straight lines are drawn to bisect the midpoint between nearest neighbor 
electrons, thereby highlighting the approximately hexagonal local symmetry. Circles are 
drawn to illustrate the approximately ring-like arrangement of electrons. The inner circles 
only pass through regular hexagons. The outer circle passes through several pentagons and 
distorted hexagons because of its proximity to the edge of the finite cluster. 

Figure 5: A particular symmetrically equivalent minimum (SEM - labelled in the text as 
fl' ) with iV near the center (c.f. Fig. 4). We are considering the limit of large N. Ring m is 
such that N >> j » 1 and it is therefore far away from the circumference of the cluster. 
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Figure 6: Typical ring m for the various spin-polarizations. Although only 12 particles are 
shown for clarity, the ring is assumed to contain a large number, i.e. N m » 1 since m » 1. 
(a) N+:N- = 1:1. (b) N+ : iV_ = 2 : 1. (c) JV+ : JV_ = 3 : 1. 

Figure 7: Theoretical lower estimate (straight line) for FQHE energy gaps obtained from 
the present theory as compared to experimental results over a range of lower mobility GaAs 
heterostructure samples. (Figure adapted from the data of Boebinger et al., Refs. 27 and 
3). Experimental data: v — | (black points), v — | (white points). As discussed in the 
text, theoretical gaps are same for v — | and | to a first approximation. 

Figure 8: Theoretical upper (dashed line) and lower (solid line - same as Fig. 7) estimates 
for FQHE energy gaps obtained from the present theory as compared to experimental results 
over a range of higher mobility GaAs heterostructure samples. Data taken from Ref. 28 
of Mallett et al. (squares) and Ref. 30 of Willett et al. (solid circles with error bars). 
Experimental data contains values for fractions v — | where p — 1,2,4,5. As discussed in 
the text, theoretical gaps are independent of p to a first approximation. 
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